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Abstract. It is shown that protostrange stars (PSSs) can 
be convective and that there are two possible scenarios de- 
scribing their turbulence. Besides the local turbulence on 
the scale which is less than the mean free path of neu- 
trinos, large-scale (~ 1 km) convection also may occur 
with properties that are similar to those of convection in 
protoneutron stars (PNSs). We thus suggest that strange 
stars can also create dynamo-originated magnetic fields 
during the deleptonization episode soon after a supernova 
explosion. Further detailed investigations are needed to 
see whether or not strange stars and neutron stars can be 
distinguished according to the differences in dynamo ac- 
tions in strange quark matter and in neutron matter. The 
magnetic fields of strange stars and neutron stars may also 
behave very differently during the accretion-phase when 
the fields decay. 

Key words: pulsars: general — stars: neutron — elemen- 
tary particles — supernovae: general — magnetic fields 



1. Introduction 

As ordinary nucleons (neutrons and protons), which are 
composed of quarks, are squeezed tightly enough at high 
temperature and pressure, they may turn into a soup of 
deconfined free quarks. This type of phase transition might 
be expected in the center of a neutron star which consists 
mainly of neutron matter. In the conventional picture of 
neutron stars with quark matter, the neutron (hadron) 
and quark phases are separated by a sharp boundary with 
density discontinuity; such stars are called "hybrid" stars. 
But, as first considered by Glendenning (1992) for "com- 
plex" systems, it is possible that bulk quark and nuclear 
matter could coexist over macroscopical distances in neu- 
tron stars; such stars are called "mixed" stars [see, e.g., 
Hciselberg et al. (2000) for a recent review about hybrid 
and mixed stars] . A further radical view is that the whole 
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neutron star should be phase-converted to be a strange 
star, which consists of nearly equal numbers of u, d, and 
s quarks and associated electrons for charge neutrality, 
if strange quark matter (SQM) is the true ground state 
of strong interacting matter (Bodmer 1971, Witten 1984, 
see also, e.g., Madsen (1999), for a recent review of the 
physics and astrophysics of strange matter and strange 
stars). Strange stars provide a sharp contrast to neutron 
stars (e.g. the mass-radius relations and the surface condi- 
tions), although both neutron and strange stars may have 
similar radii (~ 10 km) and masses (~ IAMq). 

Pulsars can be modeled as neutron stars or strange 
stars. No strong observational evidence known hitherto 
favors either of the models. It is thus of great impor- 
tance and interest in current astrophysics to distinguish 
neutron stars and strange stars observationally. Recently, 
at least three hopeful ways have been proposed for iden- 
tifying a strange star: 1, hot strange stars (or neutron 
stars containing significant quantities of strange matter) 
may rotate more rapidly since the higher bulk viscosity 
(Wang & Lu 1984, Madsen 1992) of strange matter can 
effectively damp away the r— mode instability (Madsen 
1998); 2, the approximate mass-radius (M — R) relations 
of strange stars (M oc i? 3 ) are in striking contrast to those 
of neutron stars [M oc R~ 3 ). Comparisons of observation- 
determined relations with theoretical relations may thus 
determine whether an object is a neutron star or a strange 
star (Li et al. 1999); 3, It is possible to distinguish "bare" 
polar cap strange stars from neutron stars via pulsar mag- 
netospheric and polar radiation because of the striking dif- 
ferences between the polar surface conditions of the two 
types of stars (e.g., clear drifting pattern from "antipul- 
sars", Xu, Qiao & Zhang 1999, Xu, Zhang & Qiao 2001). 
In this paper, we shall study possible differences in the 
processes of magnetic field generation by dynamo action 
in strange stars and in neutron stars which could even- 
tually lead to observational distinctions between the two 
types of stars. 

Magnetic fields play a key role in pulsar life. Unfor- 
tunately, there is still no consensus on the physical ori- 
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gin of strong fields. Naively, it is supposed that pulsars' 
fields are the result of the conservation of magnetic flux 
during supernova collapses. This idea faces at least two 
problems (Thompson & Duncan 1993, here after TD93): 
1, dynamo action can occur in the convective episodes of 
both the iron core and the newborn star; 2, very strong 
fields (> f 12 G) of pulsars are probably not "fossil fields" 
since only a few percent of white dwarfs have fields in ex- 
cess of 1 6 G. Based on the Newtonian scalar virial theo- 
rem, one can estimate the limiting interior magnetic field 
B max ~ 10 18 MiR 6 2 G of a star with mass M = Mi M© 
and radius R — R 6 fO 6 cm (Lai & Shapiro 1991). Both 
the estimate based on flux conservation and that based on 
the virial argument, do not depend on the detailed fluid 
properties of the stellar interiors. Thompson & Duncan 
(TD93, Duncan & Thompson 1992) have extensively con- 
sidered turbulent dynamo amplification in protoneutron 
stars (PNSs) and in the progenitor stars, and found that 
multipolar structure fields as strong as 10 16 G ("magne- 
tar") can be generated by PNS convection. This sugges- 
tion seems to have been confirmed by recent observations 
(e.g., Kouveliotou et al. 1998). Furthermore, there arises 
a possible way to distinguish strange stars from neutron 
stars since the dynamo-originated fields depend on the 
detailed fluid properties of the stellar interiors. One in- 
teresting question thus is: Can strange stars create mag- 
netic fields as strong as 10 16 G? (i.e., can strange stars 
act as magnetars?). If strange stars can, is there any dif- 
ference between the processes of field generation in PNSs 
and that in "proto-" strange stars (PSSs)? Theoretically, 
there is little work known hitherto on the generation of the 
magnetic fields of strange stars. It is supposed that spon- 
taneous magnetization could result in the generation of 
a compact quark star (Tatsumi 2000). Alternatively, we 
suggest that dynamo-amplification of the magnetic field 
could also play an important role in PSSs with high tem- 
peratures. In this work, we are trying to find whether there 
are any differences between the dynamo-originated fields 
of neutron stars and that of strange stars, in order to con- 
tribute to the debate on the existence of strange stars. In 
our discussion, we presume that the strong magnetic fields 
of strange stars originate in the strange cores, rather than 
in the crusts with mass ~ 1O _5 M0, since strange stars 
produced during supernova explosions cannot have such 
crusts (Xu et al. 2001). 

One difficult issue in the study of dynamo-originated 
fields of strange stars is to determine whether color su- 
perconductivity (CSC) occurs in PSS since the dynamo 
mechanism may not work effectively in a superconduct- 
ing plasma^. Recent calculations, based on a model where 
quarks interact via a point-like four-fermion interaction, 
showed that the energy gap A of zero-temperature strange 



1 In case of CSC, the field creation and evolution may also be 
influenced significantly if there exist fast dynamos with mag- 
netic Reynolds number IZm — * oo in PSSs. 



matter could be 10-100 MeV for plausible values of the 
coupling (Afford et al. 1999); thus the critical temperature 
of forming possible quark Cooper pairs T c ^A/2^5 — 50 
MeV. However, T c could be altered if a more realistic 
quark-quark interaction is used and/or if the the trapping 
of neutrinos in PSSs is included. It is also hard to suffi- 
ciently determine the temperature T(PSS) of PSSs. We 
therefore simply assume in this paper that T(PSS) > T c 
in the first few seconds (the time-scale for dynamo action 
of PSS) after PSS formation. 

2. The birth of strange stars: differential rotation 
&: convection? 

The time-scale of PNSs (or PSSs) is of the order of a few 
seconds, which is three orders of magnitude longer than 
the dynamical one. It is thus worth studying the dynami- 
cal bulk evolution in these stars, such as differential rota- 
tion and convection. Unfortunately, no numerical model of 
supernova explosion known hitherto has included the con- 
version from PNSs to PSSs, although such a conversion 
may help in modelling the burst process [e.g., 1, helping 
to solve the present energetic difficulties in getting type 
II supernova explosions (Benvenuto & Horvath 1989); 2, 
giving a reasonable explanation for the second peak of neu- 
trino emission in SN 1987A (Benvenuto & Lugones 1999)]. 
The absence of a numerical model is caused by the lack of 
(1) a full theory determining the conditions at which the 
quark matter phase transition occurs and of (2) a detailed 
understanding of the complex burning process of neutron 
matter into strange matter. 

Nevertheless, some efforts have been made in trying 
to understand the transition and combustion processes in 
detail. The phase conversion occurs in two steps: first neu- 
tron matter deconfinement occurs on a strong interaction 
time scale ~ 10 -23 seconds, then chemical equilibration 
of the deconfined quark matter takes place on a weak in- 
teraction time scale ~ 10 -8 seconds. Additional neutrinos 
and energy are produced in the second step (Dai, Peng & 
Lu 1995). Further calculations of such transitions (Anand 
et al. 1997; Lugones & Benvenuto 1998) have also consid- 
ered the effect of strong interactions, the effect of finite 
temperature and strange quark mass, and the effect of 
trapped neutrinos. Recently, Benvenuto & Lugones (1999) 
explored the occurrence of deconfinement transition in a 
PNS modeled by Keil & Janka (1995) and found that the 
deconfinement appears as long as the bag constant B < 
126 MeV fm -3 . Various estimates of the bag constant in- 
dicate that the preferred value of B lies in the range of 
60MeV far 3 < B < HOMeV far 3 (Drago 1999), which 
means that deconfinement is very likely to happen. Such 
2-flavor quark matter may be transformed immediately 
into a 3-flavor one if SQM is absolutely stable. From a 
kinetic point of view, Olinto (1987) has calculated for 
the first time the conversion of neutron stars into strange 
stars, suggesting a deflagration mode with a burning ve- 
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locity range from 10 4 km/s to a few cm/s. However it is 
found (Horvath & Benvenuto 1988; Benvenuto, Horvath 
& Vucetich 1989; Benvenuto & Horvath 1989) that such 
slow modes are unstable. This instability would be self- 
accelerated, and the burning should occur finally in deto- 
nation modes, although the transition from deflagration to 
detonation has not been well understood (Lugones, Ben- 
venuto & Vucetich 1994). In order for the combustion to 
be exothermic, there exists a minimum density p c for det- 
onation to be possible, p c ~ 2 times of nuclear density for 
some cases (Benvenuto & Horvath 1989). Thus the deto- 
nation flame can not reach the edge of the compact core, 
and the outer part of the PNS would be expelled. As a re- 
sult, a strange star is formed with an almost "bare" quark 
surface which is essential to solve completely the "binding 
energy" problem in some current pulsar emission models 
(Xu ct al. 1999, Xu et al. 2001). Both the extra neutrino 
emissivity and the detonation wave can favour a successful 
core-collapse supernova explosion. 

2.1. Differential rotation 

It is expected that PNSs (thus PSSs) have a strong dif- 
ferential rotation (e.g., Janka & Mdnchmcycr 1989, Gous- 
sard et al. 1998). This issue is uncertain, however, due 
to the lack of a rotating core model in the pre-supernova 
evolution simulations. The iron core collapses, being trig- 
gered by electron capture (for lower entropy core) and/or 
photodisintegration (for higher entropy core), almost in 
the same way as for free fall. Many initial core models 
of highly evolved massive stars without rotation have ap- 
peared in the literature (e.g., Arnett 1977, Bruenn 1985), 
in which the density po - radius ro relation can be fitted 

by 



A)Oo) 



i=l 
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where poi, a\ are constants. For the first order approxima- 
tion, we let n = 1 in our following discussion. In this case, 
the Lagrangian mass coordinate M(rn) is 



M(r ) 



Airp 



l[2-(2 + 2a 1 r + a^)e- airo ] 



(2) 



Setting A"f(10 8 ) - 1.4M Q and p Q {W 
for the core with radius Ro = 10 8 



10 s g cm" d 
one obtainsp| 

poi = 2 x 10 io g cm" 3 , cei = 5.4 x 10" 8 cm" 1 , by Eq.([y) 
and Eq. (0) . We assume that the core before collapse with 
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2 The density model p (r ) = p ie~ air ° with p 01 = 2 x 10 10 
g cm -3 , ai — 5.4x 10~ 8 cm -1 , based on which we will calculate 
the nature of differential rotation of PSSs (or PNSs), represents 
well the initial core model of Bruenn (1985, Fig.l, there 7?o ~ 
10 8 cm, X(10 8 ) ~ 1.4M , p (0) = 2 x 10 10 g cm" 3 , and 
po(10 8 ) = 8 x 10 7 g cm" 3 ). We thus think that the computation 
according to this simplified density model can give the correct 
order of magnitude result for the actual collapsed rotating core. 



a uniform rotation of period P can a l so be approximated 
by Eq.(l) for lack of a rotation core model in simulations, 
and we investigate below the collapse of the iron core with 
the above parameters. Here we also assume that pulsar 
rotation is mainly the result of angular momentum con- 
servation during the core collapse process, rather than of 
the kick at birth (Spruit & Phinney 1998)0. 

Based on the approximation, the total rotational en- 
ergy Eq of the core is 



En = 



167T 3 p i 3(ai,i?o) 



3P 2 



and the total angular momentum Mq is 
167T 2 p i 3(ai,i? ) 



M 



3P 



(3) 



(4) 



where the function 3(a, R) is defined by 

S(a,R) = 24-(24+24aP+12a 2 P 2 +4a 3 P 3 +a 4 P 4 )e" afl (5) 

It is found that lim a ^o^s(a,R)/a 5 = i? 5 /5. Contrary to 
the case of the core before collapse, newborn strange stars 
may be well approximated as objects with homogeneous 
density (Alcock et al. 1986). For a strange star with mass 
M = 1.4M© and radius R = 10 6 cm, we find the den- 
sity p ~ 7 x 10 14 g cm" 3 . During the adiabatic collapse 
process, in which the entropy of each mass element does 
not change significantly, toroidal forces are negligible (al- 
though a poloidal force can cause a shock wave). The an- 
gular momentum of each mass element is thus conserved, 
and PSSs should be in differential rotation. We assume 
that the material at a shell with radius ro of the iron core 
before collapse contracts to another shell with radius r in 
a PSS, 



r 3 = 4^M(ro). 



(6) 



According to angular momentum conservation of each 
mass element, one arrives at the velocity field V dependent 
on r and 9 (polar angle), 



V = 



27rsin# 



1° 
r 



(7) 



where ro is a function of r through Eq. (^) , and the mass 
conservation law Airr^podro = AiTr 2 p&r has been used. 
Owing to the momentum transport by neutrinos, magnetic 
fields, and turbulence, a uniform rotation is approached 
after a certain time and the final rotational energy Ef can 
be estimated as 



Ek = ^PR 5 ~ 2.3 X 10 46 P" 2 , 
corresponding to the angular momentum Mf, 
16tt 2 



M f = 



15P 



pR* 



(8) 



(9) 



3 The recent observation of spin-kick alignment of Crab and 
Vela pulsars supports this suggestion. 



4 



R.X. Xu & F.H. Busse: The birth of strange stars and their dynamo-originated magnetic fields 



where P is the period of newborn uniformly rotating 
strange star. According to angular momentum conserva- 
tion in the collapse process, M = Mf , the rotation period 
Pq of the initial iron core before collapse is 



Po 



5/9 iQ ! (ai, Rq) 



P ~ 4.68 x 10 P, 



(10) 



for typical parameters, i.e., the core rotates with a period 
of about 50 seconds for an initial pulsar period P ~ 10 ms. 
From Eq.(0), one obtains the velocity derivative |VV|, 



(11) 



as a function of 9 and r. 

Fig.|l| shows the velocity and the velocity derivative of 
this differential rotation scenario for P = 10 ms, respec- 
tively, based on Eq.(0) and Eq.@. We find |W| ~ 10 4 
s~ 4 in the outer part of PNSs or PSSs. Such differential 
rotation may play an important role in the creation of 
magnetic fields of PSS or PNS, as will be discussed in 
section 3. 

Since the dynamo actions are in the outer convective 
layer with thickness L ~ l p ~ 10 5 cm (see Eq.(|l2]) and 
section 3), we give the velocity difference of differential 
rotation, AV = V(r = 10 6 cm) - V(r = 9x 10 5 cm), as a 
function of P in Fig.||. 

2.2. Convection 

Turbulent convection in PNSs has been extensively inves- 
tigated before (e.g. Burrows & Lattimer 1988, Wilson & 
Mayle 1988, Miralles, Pons & Urpin 2000). In the Kelvin- 
Helmholtz cooling phase, both the negative gradient of the 
entropy and of the lepton fraction can drive convection in 
newborn neutron stars, given that the mean free paths of 
leptons are much smaller than the convection length scale. 
Protostrange stars form after the SQM phase-transition 
with a time scale of ~ R/c ~ 10 -4 s, which is much shorter 
than that of neutrino diffusion and thermal evolution, as- 
suming that the actual combustion mode is detonation 
(Benvenuto, Horvath & Vucetich 1989). 

We expect that PSSs are convective, since negative 
gradients of entropy and neutrino fractions would also ap- 
pear in the outermost layers which can lose entropy and 
neutrinos faster than the inner part. Following arguments 
used by Thompson & Duncan (TD93) in the case of PNSs, 
we compare the radiative and adiabatic temperature pro- 
files to see whether convection occurs in PSSs. We assume 
neutrinos in the outer part of PSS are nondegenerate since 
there the neutrinos are lost rapidly. Iwamoto (1982) has 
shown that for the mean free path I of the scattering of 
nondegenerate neutrinos by relativistic degenerate quark 
matter with temperature T, l/l cx T 3 holds. Thus, the 
neutrino opacity in SQM scales as T 3 (rather than T 2 as 
for the case of PNSs), and the radiative temperature pro- 
file in PSSs is much steeper, T(r) <x p(r). We can employ 



the analytical equation of state for SQM in the case of zero 
strange quark mass (m s — 0) and zero coupling constant 
(a s = 0) of strong interaction (Cleymans et al. 1986; Ben- 
venuto et al. 1989) to estimate the adiabatic temperature- 
pressure relation. The entropy per baryon is S ~ 7>n 2 T j p 
for the case of chemical potential p ~ 300 MeV and tem- 
perature T ~ 30 MeV. Therefore T(r) ~ (p(r) + B) 1/4 
for SQM moving with fixed S is much less steep, and 
we can thus expect that a negative entropy gradient ap- 
pears in the outer layer which is unstable to convection, 
similar to the case of PNSs. Accordingly we suggest that 
Schwarzchild convection exists in PSSs, whereas it will de- 
pend on detailed simulations whether Ledoux convection 
in PSSs can be established. If the timescale tnm for neu- 
tron matter convection is much smaller than the neutrino 
diffusion time (a few seconds) , we may expect that Ledoux 
convection takes place in a PSS since the negative gradi- 
ent of the lepton fraction is nearly the same for PSSs as 
that for PNSs. 

In the following, we try to estimate the properties of 
convection in PSSs. The local pressure scale height l p in 
PSS is 



pg 



'' 2 x 10 5 cm, 



(12) 



wher 



e g 



GM 



= 1.33 x 10 14 Mi i?g 2 cm s~ 2 (M x 



M/Mq, i?6 is the radius R in 10 6 cm) is the gravitational 
acceleration, p = (pc 2 — 4B)/3 ~ 10 34 dyne cm~ 2 is the 
typical pressure in the outer convective layer, and the typ- 
ical density there being chosen as p — 5 x 10 14 g cm~ 3 
(Alcock et al. 1986). We can assume the thickness L of 
the outer layer where dynamo action exists to be equal to 



this scale height, L 
VT thus is 



l p , and the temperature gradient 



VT 



T 
L 



5 x 10 5 Tn K/cm, 



(13) 



where Tn is the SQM temperature in 10 11 K. The numer- 
ical equation of state of SQM including the effect of non- 
zero strange quark mass [Eq.(29) of Chamaj & Slominski 
1989] is 

p = 1.07 x 10 14 i?60 + 3F/c 2 + 1.02 x i(r 9 P - 6455 

+ 1.34 x i -i4r 2 P 09619 , (14j 

where P/c 2 = p/c 2 + 1.07 x 10 14 fi 60 (B 6 o is the MIT 
bag constant in 60 MeV fm -3 , p is the external pressure). 
According to Eq. ( |l4"| ) the coefficient of thermal expansion 
a is 

a = --(^) = 2.68 x lO-^TuP 09619 ^ 1 K" 1 . (15) 

a = 1.33 x 10~ 14 KV 1 for T = 10 n K and p = (Usov 
1998). 

There are two factors which can cause viscous stresses 
in a PSS: neutrino transport and quark scattering. 
Neutrino-induced viscosity dominates in a PSS on scales 
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Fig. 1. The calculated velocity and velocity deriva- 
tive profiles for different polar angles 9. Dot- 
dashed, dotted, dashed, and solid lines are for 9 — 
7r/8, 7r/4, 37r/8, 7r/2, respectively. The lower (upper) 
4 lines, scaled by the left (right) ordinate, are for 
velocity (velocity derivative). P = 10 ms in the cal- 
culations. 



Fig. 2. The velocity difference in the convective 
outer layer with thickness of L ~ 1 km as a function 
of pulsar initial period P for different polar angles 9. 
See Fig]j] for the definition of lines. 



that are large compared to the neutrino mean-free path 
We use the neutrino mean- free path of nondegenerate 
neutrino scattering in SQM for a dimensional estimate 
(Iwamoto 1982), 



1 = 1.7 x 10 2 pJ /3 E- 



3 

100 



cm, 



(16) 



where pis is SQM density in 10 15 g cm~ 3 and Exoo is 
neutrino energy in 100 MeV. According to Eq.(ll) of 
TD93, the neutrino mean free path in nuclear matter, 
Z N - 10 2 pis /3 cm for T = 30 MeV, which is of the same 
order as I. Unfortunately, no well-determined neutrino vis- 
cosity in PSS has appeared in the literature. We thus just 
estimate the neutrino- induced viscosity by a simple kinetic 
argument (Wilson & Mayle 1988) 



v = \lc£, = 1.7 x 10 10 p 15 2/3 £r 00 



1 — 1 

cm s . 



(17) 



where we have assumed the ratio £ of the neutrino en- 
ergy density to the quark one to be ~ 10~ 2 . However, 
the kinematic viscosity due to quark scattering in SQM 
(Heiselberg & Pethick 1993) is much smaller 



0.1 pT K/O.I)- 5 / 3 T n 5/3 



cm 2 s . 



(18) 



Therefore turbulent convection may have a scale I, 
which finally should be damped on a small scale by v q . 

Two possibilities arise for the scenario of turbulence in 
PSSs. The first one is that the neutrino viscosity can ef- 
fectively inhibit a large-scale convection with length scale 
L, and local convection with scale I will exist. In this 
case the local thermal diffusivity due to quark scattering 
in SQM (Heiselberg & Pethick 1993) is 



Kq ~0.39 p% 3 (a s /0.1) 
The Prandtl number is 



P iq = ~ 0.25, 



the Rayleigh number is 



11 



2 — 1 
cm s . 



Raq — 



agWTl 4 _ n n ^ 16 



7.2 x 10 1 



KqVq 

and the Coriolis number is 

2m 2 o 
Tq = 2 x 10 s . 



(19) 



(20) 



(21) 



(22) 



6 



R.X. Xu & F.H. Busse: The birth of strange stars and their dynamo-originated magnetic fields 



Malkus (1954) had estimated the mean-square values (v 2 ) 
of the fluctuating velocity from the Boussinesq form of the 
hydrodynamic equations, and found [Eq.(64) in Malkus 
(1954)] 



K (U 



Pa 



Nl/2 



(23) 



where k is the thermal diffusivity, d is the length scale 
of convection, P a is the Rayleigh number, and R ac is the 
critical Rayleigh number. This relation was supported by 
experiments [Fig.2 in Malkus (1954)] for P a = 10 5 - 10 9 . 
Also, Clever & Busse (1981) found the perturbation en- 
ergy (v 2 ) increases nearly proportional to P a — R ac at 
values of R a > 10 3 for low-Prandtl-number convection 
[P r > 0.001, see Fig.ll in Clever & Busse (1981)]. Accord- 
ing to Eq.(|23|), the turbulent convective velocity on this 
scale is 



(24) 



R 



at) ■ 



v q ~ 3.5 x 10 5 cm/s, 

if we choose k — K q , d = I, R a 

The second scenario is that the neutrino viscosity is 
not high enough to inhibit the large-scale convection, and 
convection with scale L is possible. We can also obtain 
some dimensionless numbers for this case. Since both the 
thermal energy and momentum in PSSs are transported 
by neutrinos now, the thermal diffusivity k can be esti- 
mated to be (Wilson & Mayle 1988) 

k = \lc = 1.7 x 10 12 pi5 2/3 ^roo cm 2 s-\ (25) 
by a simple kinetic argument. Thus the Prandtl number 



P t = - = £ ~ 0.1 - 0.01. 



(26) 



For a rotating spherical fluid shell with thickness L and 
angular velocity f2, the Rayleigh number is 



R a — 



agVT£4 ~ 6.0 x 10- 23 T L 3 



(27) 



i.e. R a ~ 6.0 x 10 3 for T ~ 10 11 K, p = and L ~ 10 5 
cm. The Taylor number r 2 can be calculated by 



2fl L 2 



= 1.2 x 10~ 10 p% 3 Ef 00 n L 2 . 



(28) 



r - 1.2 x 10 3 for Q ~ 10 3 s^ 1 and L ~ 10 5 cm. The critical 
Rayleigh number R ac can be estimated (Zhang 1995) by 



Pa 



10 4 



(29) 



Turbulent convection is thus possible for R a > R ac with a 
velocity which can be estimated to be (d = L) 



v & 4 x 10 8 cm/s, 



(30) 



according to Eq.([23|). Alternatively if the mixing-length 
prescription (Bohm-Vitense 1958, TD92) for a fluid of 



semidegenerate fermions is adapted, the convective veloc- 
ity in PSS becomes 



1 L 



"ml 



r >2 



2r 4irR 2 p 



) 1 / 3 ^ 6 .8xl0 7 L^p-5 



1/3 



where V = SSE = 

dlnp pc 2 -4B 



PC 



5, L52 is the convective 



luminosity in 10 52 erg s 1 . We note that v m \ and v are 
nearly of the same order, while for PNSs with T = 5/3, 
the convective velocity is 5.4 x 10 7 L X J^ p 15 1 ^ 3 cm s _1 ~ 10 s 
cm s^ 1 for p ~ 5 x 10 14 g cm -3 . Thus both, the convection 
in PSSs and in PNSs have convective velocities of about 
10 s cm s _1 if large-scale convection is possible. Also the 
overturn timescales of PSSs and of PNSs are of the same 
order, r con = l p /v ~ a few ms, and the Rossby number is 
R Q = P/r con ~ 10 3 Pi in this possibility (Pi is the value 
of the initial period in seconds). 

Another cause of turbulent motion in PSS (and PNS) 
is differential rotation. The Reynolds number 1Z of this 
shear flow is 



11 



\7UL 2 



(32) 



where Vt7 and L are typical velocity derivative and length 
scales, respectively. For the local convection with scale ^ /, 



n 



10 . For the convection with scale L, 1Z ~ 5.8 x 



10 . According to the experiment by van Atta (1966, Fig. 3 
there), turbulence requires Reynolds number 1Z > 10 4 , 
thus in our case small scale (^ I) turbulence may appear 
owing to the differential rotation. 



3. Dynamo actions in proto-strange stars 

Magnetic fields in a static plasma of finite electric con- 
ductivity o are subject to diffusion and dissipation, the 
timescale of which is r i im = L 2 / K m — 4tt<jL 2 /c 2 (ft m = 
c 2 /(47rcr) is the magnetic diffusivity). However, in a flow- 
ing plasma, magnetic fields may be created with certain 
velocity fields. This is the so-called "dynamo" process, 
which converts kinetic energy into magnetic energy (e.g., 
Moffatt 1978). The equation describing the time- variation 
of the magnetic field B governs the dynamo action 



SB 

~dt 



V x (vxB)+K-¥B, 



(33) 



where 1Z m — T Q \ im /T v = UL/ n m is the magnetic Reynolds 
number (t v = L/U is the advection timescale). The mag- 
netic diffusivity can be estimated 



Km = -i— = 5.8X10-V 879 {a s /Q.lf/' A T'^ cm 2 s~ l (3A) 

for SQM [here a is calculated from Eq.(39) and Eq.(28) of 
Heiselberg & Pethick (1993)]. Thus K m - 10 16 for PSSs. 
But for PNSs, the magnetic diffusivity is (TD93) 



N5/3 T 5/3 a 



k£ = 1.5 x lO" 5 ^ 173 cm 2 s" 1 , 



(35) 
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for neutron matter with an electron fraction of 0.2, which 
is about two orders smaller than K m . From Eq.(0) and 
Eq.([34j), one obtains the magnetic Prandtl number P m 

P m = — = 3.0 x 10 12 pH 9 E w l (a s /0.ir 5 / 3 T n 5/3 . (36) 

Km 

3.1. An estimate of field strength 

Thompson & Duncan (TD93) concluded that the domi- 
nant kinetic energy to be converted into magnetic energy 
in the dynamo action of PNSs is the convective energy 
since the initial pulsar rotation period is probably much 
larger than the overturn time (~ 1 ms) of a convective 
cell, although large-scale a — f2 dynamo action is essen- 
tial for neutron stars with R Q ^ 1 (i.e., initial period ^ 1 
ms) to produce very high fields ( "magnetars" , Duncan & 
Thompson 1992). However, here we suggest that rotation 
can not be neglected even for pulsars with typical initial 
periods, because the differential rotation energy density 
Ed jfi El ~8x 10 31 erg cm -3 is even larger than the turbu- 
lent energy density pv 2 /2 ~ 5 x 10 30 erg cm~ 3 for P = 10 
ms. If the large-scale convection scenario discussed in sec- 
tion 2.2 is possible, most of the differential rotation energy 
may be converted to magnetic energy by dynamo action. 
Let's first estimate the total differential rotation energy 
Pd in the case when the angular momentum of each mass 
element is conserved in the collapse, 

E * = C rodr ~ 3.1 x 10 46 Pr 2 erg. (37) 

Jo 

For a typical initial period P = 10 ms, one gets ~ 10 50 
erg. The dominant energy in the core-collapse supernovae 
is the gravitational energy E g ~ 0.6 GM 2 /R ~3x 10 53 
erg for M = 1.4 M Q and R = 10 6 cm, if a strange star 
is residual. We thus find that a significant part (~ 0.03% 
for P = 10 ms) of the gravitational energy has to be con- 
verted to the rotation energy if angular momentum is con- 
served. If most of this differential rotation energy would 
be converted to magnetic energy by a — f2 dynamo action, 
ItT ~ jg> ■ one obtains the saturation magnetic field in 
the interior of the stars, 

P sat ~ 4.3 x 10 14 Pf 1 G, (38) 

which is near the value estimated by Thompson & Duncan 
(TD93) who assumed that most of the convective energy 
is converted to the magnetic one. i? sat ~ 10 16 G for a 
typical period P = 10 ms. Certainly, the actual "dipole" 
magnetic field B p of pulsars should be only a fraction of 
P sa t- For example, in the upper convection zone of the 
sun, the rms value the magnetic field is only 10~ 2 P sa t. If 
the same is assumed for the pulsars, the rms value of fields 
in PSSs or PNSs should be - 10 12 Pf 1 G. However, the 
poloidal magnetic field should be another fraction of this 
rms value, thus 

B p < 10 12 P- 1 G. (39) 



Observations for ordinary pulsars show that the fields are 
distributed from - 10 11 G to 5.5 x 10 13 G (PSR J1814- 
1744). Therefore, it is likely the initial period P of ordinary 
pulsars could be in the range of a few hundreds to decades 
of milliseconds. If P < 10 ms, much stronger dipole fields 
(e.g., ~ 10 14 — 10 15 G for magnetars) can be generated. 

In the other scenario, where diffusivities based on neu- 
trino scattering inhibit large-scale convection, the differ- 
ential rotation is damped due to the momentum transport 
by neutrinos. In this case the dynamo may be driven domi- 
nantly by turbulent convection, and the field strength can 

be estimated as B = ^Jinpv 2 ~ 10 13 G, assuming an 

equilibrium state between kinetic and magnetic energy. 
The efficiency of converting differential rotation energy 
into magnetic energy by dynamo action in this scenario 
may be much smaller than in the case where convection 
with the large scale L exists since a significant fraction of 
differential rotation energy is likely to be converted into 
thermal energy due to the high neutrino viscosity. 

Thus, dynamo action may amplify significantly the 
field before PSSs cools down to temperatures T smaller 
than the critical one, T c . When T < T c , CSC appears, and 
the field would exist as a fossil one for a very long time 
since r i im — ► oo if a — > oo. In fact, Afford et al. (2000) 
investigated recently the effect of CSC on the magnetic 
fields, and found that, unlike the conventional supercon- 
ductors where weak magnetic fields are expelled by the 
Meissner effect, color superconductors can be penetrated 
by external magnetic fields and such fields can exist stably 
on a timescale longer than the cosmic age. 

Eq.(|39"|) may have an observational consequence for 
pulsar statistics, which in turn could test Eq.(|39"|). In the 
magnetic dipole model of pulsars, (3.2 x 10 19 ) 2 P^ = B 2 
(note: fields are assumed not to decay in SQM), and the 
rotation period P is a function of time t due to energy 
loss. P(0) denotes the initial period of pulsars. Consider- 
ing Eq.(f§), one has P(t)^p- = <P(0)~ 2 , the solution of 
which is 

P{t) 2 = P(0) 2 + 2C P(0)- 2 t, (40) 

where £ ~ 4.3 x 10~ 18 s 3 . Observational tests of this equa- 
tion are highly desirable. 

The condition E& < E g gives a limit for the initial 
pulsar period P. Based on Eq.(|37]), one gets P > 0.2 ms. 
Actually, P ^> 0.2 ms since the efficiency of converting 
gravitation energy into differential rotation energy during 
collapse may be very small. It is thus doubtful that super- 
novae can produce pulsars with submilliseconds periods. 
For recycled millisecond pulsars the above estimate is not 
relevant. 

3.2. Fast dynamos 

The magnetic field amplification processes in newborn 
pulsars are essentially fast dynamos because of the high 
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magnetic Reynolds numbers (TZ m ~ 10 16 — 10 17 for large- 
scale convection, TZm ~ 10 10 for local turbulence) of both 
PSSs and PNSs. Unfortunately, the question of whether 
or not fast dynamo exist has not been answered theo- 
retically although many numerical and analytical calcu- 
lations strongly support the existence of kinematic fast 
dynamos for given sufficiently complicated flows (Soward 
1994, Childress & Gilbert 1995). It is worth studying fluids 
without magnetic diffusion since the diffusion timescale is 
much longer than the advection timescale, r i im S> r v , 
for fast dynamos. The complex flows, such as stretch- 
twist-fold, may effectively amplify the field in this case. 
The amplified strong magnetic fields are concentrated in 
filaments with radii if, which can be estimated to be 
k ~ TZ m 1/2 L ~ 0.01/ % /7^~u7 mm {TZ m = TZ nll6 x 10 16 ) 
by equating the diffusion timescale of filament field, If / K m , 
to the advection timescale r v . Fast dynamos in PNSs have 
been considered by Thompson & Duncan (TD93) who sug- 
gested PNS dynamos as the origins of pulsar magnetism. 
As the fluid parameters (e.g., 7?. m , v, k) of both PSSs and 
PNSs are similar, fast dynamos may also work for new- 
born strange stars. 

There are three timescales in the fast dynamo of 
PSSs: the diffuse timescale r h m = L 2 /n m , the advection 
timescale r v = L/U, and the buoyancy timescale Tb, 

Tb * Ac \[^M RLl/2pl/2B ~ 1 ~ lOnisiW^X^l) 

(R = R# x 10 6 cm, L = L h x 10 5 cm, B = B 16 x 10 16 
G). Bi e ~ 1.12 p¥ 5 2 v$ (v = Vg x 10 8 cm s _1 ) if we assume 
that kinetic and magnetic energy densities are in equilib- 
rium. Filaments can be generated inside the convective 
layer since turbulent convection is violent on small scales 
(~ 100 cm). t v ^ Tb <C T hm, which means that filaments 
rise to the stellar surface by magnetic buoyancy as soon 
as strong fields are created by fast dynamos. Magnetic 
buoyancy flow or convection could create and amplify the 
poloidal (transverse) field from the azimuthal (aligned) 
field, at the same time, the aligned field could also be cre- 
ated and amplified from the transverse field by differential 
rotation as discussed in section 2.1. Thus these processes 
may complete a "dynamo cycle" in newborn pulsars. 

Let's give some estimates for fast dynamos. There are 
mainly two types of flows in PSSs: eddy (convection) and 
shear (differential rotation) . Actually, these flows are cou- 
pled. However, we may deal with them separately in order 
to have an overview of the field generation processes. For 
pure straining motion with velocity field V, 

V = r-\-x,y,0), (42) 

which represents eddy flow to some extent, the magnetic 

1 1% 

field can be amplified considerably to -B£ lax = BolZrA at a 
timescale i^ ax = jT v \nR,]l 2 (Soward 1994). For an initial 
field strength B ~ 10 10 G, B^ ~ 10 18 ^16 G and 



*max ~ 18r V 3ln7^^ 6 ms (r v = r v3 x 10 3 s). On another 
hand, the linear shear flow (Soward 1994), 

V = T v - 1 ( 2 /,0,0), (43) 

can amplify the field to B s mlix ~ B Q TlU 3 - 2 x lO 15 ^ 7 ^ 
G at a growth time tf nax ~ T V TZ][ 3 ~ 2\T V ^R^ & s for 
B = 10 10 G (T v = 10- 4 T v4 ~ L/(V(R)-V(R-L)), 
see Fig.|J). The timescale for field generation is thus from 
10" 5 s (for 7l m ~ 10 10 ) to 10 s (for K m ~ 10 16 ), while the 
amplified fields could be 10 12 (for K m ~ 10 10 ) to 10 18 (for 
7Z m ~ 10 16 ) G in filaments. The magnetic fields emerging 
from under the stellar surface are likely to be much smaller 
than -B m ax in filaments because the filaments will expand 
to an approximately homogeneous field near and above the 
surface owing to magnetic pressure. Incoherent or coherent 
superposition of this buoyant field flux might result in the 
observed dipole moment of pulsars (TD93). 

Differential rotation may play an essential role for the 
generation of large-scale magnetic fields through a — Q 
dynamo process, but even in the absence of differential 
rotation, large-scale magnetic fields may be created. Ow- 
ing to the alignment of small-scale convection rolls parallel 
to the axis of rotation, global magnetic fields can be gen- 
erated as has been shown in various dynamo models (see, 
for instance, Busse 1975). 

4. Any differences of dynamo-originated fields? 

As discussed in the previous section, magnetic dipole fields 
stronger than 10 11 G could be generated by dynamos in 
PSSs. However, if pulsars are born as strange stars, an 
acute question is how to explain the fields of millisecond 
pulsars with ~ 10 8 G since fields don't decay in strange 
stars with CSC. Millisecond pulsars are supposed to be 
spun- up ( "recycled" ) as the result of accretion from com- 
panion stars in their histories. They could be strange stars 
with ~ 10 _5 Mq crusts, although a lot of ordinary pulsars 
might be strange stars with bare polar caps and much 
thinner crusts (Xu et al. 2001). We suggest that accre- 
tion may reduce the dipole field strength from ~ 10 12 
G to ~ 10 8 G of strange stars in four possible ways: 
(1) the accreted matter may screen or bury the origi- 
nal fields (Bisnovati-Kogan & Romberg 1974, Zhang et 
al. 1994, Zhang 2000); (2) field decays in the accretion- 
heated crust (Ronar & Bhattacharya 1997); (3) the large 
scale fields (e.g., dipole structure) could be changed sig- 
nificantly because the plasma accreted onto the polar caps 
would squeeze the frozen magnetic fields toward the equa- 
tor (Cheng & Zhang 1998), and the dipole field thus ap- 
pears to decay; (4) magnetic fields frozen in the bottom 
crust may dissipate and be annihilated during the com- 
bustion process of the bottom matter into SQM when the 
crust becomes too heavy to be supported by the Coulomb 
barrier. 

A notable difference between strange stars and neu- 
tron stars is that the magnetic field of a strange core is 
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stable and does not decay due to CSC, while the field 
of a neutron star should be expelled from the interior to 
the crust where the Ohmic dissipation occurs. In present 
models for neutron star magnetism, the fields permeate 
cither the whole star or only the crust (see, e.g., a short 
review by Mitra et al. 1999). It is found by observational 
and statistical analyses that a pulsar's field can only de- 
crease substantially in the accretion-phase but does not 
decay significantly during the pulsar lifetime (Hartman et 
al. 1997). New calculations of the timescale for Ohmic dis- 
sipation in the crust have shown that the field can persist 
for more than 10 10 years (Sang & Chanmugam 1987). It 
thus may be difficult to distinguish the field evolutions of 
strange stars and neutron stars in the radio pulsar-phase. 
Nevertheless, the field decay modes in the accretion-phase 
could be quite different for neutron stars and strange stars. 
In fact the items (l)-(3) in the above proposals for field 
decay are relevant to both neutron stars and strange stars, 
but item (4) can only apply to strange stars. Furthermore, 
items (l)-(3) could result in different physical processes 
for strange star and neutron star because they have very 
different structures. For example, as matter accretes, the 
radius of a strange star increases, while the radius of a 
neutron star decreases. Also, the material in neutron star 
crust moves continuously, and the movement just pushes 
and squeezes the original field. However for strange stars, 
a SQM phase-transition occurs when the accreted ma- 
terial moves across the Coulomb barrier with thickness 
of ~ 200 fm. Actually, this phase-transition process has 
been included in the explanation of some burst phenom- 
ena, such as bursting X-ray pulsar GRO J1744-28 (Cheng 
et al. 1998), but the field evolution in accreting strange 
star has not yet been discussed extensively. We propose 
here that strange stars and neutron stars may be distin- 
guished by their field decays during the accretion-phases. 
Future investigations of this issue will be of interest. 

As addressed in the previous sections, many fluid pa- 
rameters for the large-scale convection scenario in PSSs 
and PNSs are similar. Both kind of stars have convective 
layers with thickness ~ 10 5 cm and flow velocities ~ 10 8 
cm/s. The scale length of turbulent eddies of both PSS 
and PNS is ~ 2 m since neutrino viscosity damps flow 
only on scale greater than I ~ I N ~ 2 m. The neutrino 
viscosities and the magnetic Reynolds numbers in PSSs 
and PNSs are also not quite different: v ~ 10 10 cm 2 /s, 

- 10 8 cm 2 /s; K m ~ 10 16 , 11™ ~ 10 17 . Therefore, the 
general configurations and strength of dynamo-generated 
magnetic fields in PSSs and PNSs may be similar. How- 
ever, the fluid properties calculated in this paper for SQM 
are rather uncertain since our knowledge of SQM fluid is 
less than that of nuclear matter. Detailed studies of the 
fluid properties in PSSs, such as the neutrino fraction and 
viscosity, are necessary and may help us to see differences 
in the dynamo-originated fields. In addition, the timescale 
available for dynamo action in PSSs may be significantly 
smaller than that in PNSs since the energy gap for proton 



superconductivity is order of 1 MeV (rather than 10-100 
McV in SQM). This property may also affect the structure 
of the fields. 

5. Conclusion and discussion 

Pulsars could be neutron stars and/or strange stars. The 
origin of neutron star magnetic fields has been discussed in 
literature, while few papers are concerned with the origin 
of fields in strange stars. In this paper, we have inves- 
tigated the dynamo action in PSSs, and have suggested 
that strange stars can have magnetic fields of dynamo ori- 
gin similar to those in neutron stars. Our main conclusions 
are as follows. 

1. A significant fraction of gravitational energy has to 
be converted to differential rotation energy if the angular 
momentum of each mass element is conserved in the core 
collapse. Assuming an initial rotating core to be approxi- 
mated by the model of Bruenn (1985), we have calculated 
the rotation velocity, the velocity derivative, and the dif- 
ferential rotation energy of a nascent strange star. 

2. It is found that PSSs may have a convective layer 
of thickness ~ 2 km in the outer part. The large-scale 
convection has a velocity ~ 10 8 cm/s, while local turbulent 
eddies have a scale of ~ 1 m and a velocity of ~ 10 5 cm/s. 

3. The energy density of differential rotation could be 
larger than the turbulent energy density if the pulsar ini- 
tial period P & 10 ms. Assuming that most of the dif- 
ferential rotation energy is converted to magnetic energy 
by dynamo action, we obtain the dipole field strength as 
a function of P and the pulsar period evolution due to 
magnetic dipole radiation. 

4. The fields (10 12 - 10 18 G) amplified by fast dynamo 
action are concentrated in filaments with initial radii ~ 
0.01 mm - <~ 1 cm and with growth times 10 -5 s to 10 s. 

5. Strange stars and neutron stars are expected to have 
different accretion-induced field decay processes, which 
could be used to distinguish them in the future. 

6. Convection with the large scale L is less likely to 
exist in PSSs than in PNSs. 

It is currently believed that anomalous X-ray pulsars 
(AXPs) and soft gamma repeaters (SGRs) are magnetism- 
powered in the "magnetar" model (Duncan & Thompson 
1992, Thompson & Duncan 1995). The magnetic recon- 
nection near the surface causes the emission in AXP and 
SGR. Duncan & Thompson (1992, TD93) suggested that 
the key parameter that determines whether a PNS be- 
comes an ordinary pulsar or a magnetar is the Rossby 
number. If magnetars are strange stars, we propose that, 
besides this key parameter, the initial temperature and the 
density of trapped neutrinos can also affect the formation 
of a magnetar since dynamo action could not be very ef- 
fective when CSC appears. Strange stars with very strong 
magnetic field may act as "magnetars" since they could 
have similar differentiated structure as neutron stars (Bcn- 
venuto et al., 1990). However the neutron magnetar model 
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faces a crisis indicated by Perez Martinez et al.(2000) in 
that neutron stars may undergo a transverse collapse if 
their fields exceed a critical value. But for strange mag- 
netars, a relatively small magnetic momentum and large 
chemical potential of free quarks may favor the forma- 
tion of very high fields, although a further investigation 
is needed to find the critical field strength beyond which 
a strange star cannot be sustained against transverse col- 
lapse. 

This paper may have two implications for the studies 
of the r— mode instability (Anderson 1998; Madsen 1998, 
2000; Lindblom et al. 1998). (1) As discussed in section 
2.1, PSSs and PNSs should rotate strongly differentially 
(particularly for P ^ 10 ms). But the numerical studies 
that have appeared in the literature are for stars with uni- 
form rotation. Certainly, it is reasonable to assume that 
the result for uniform rotation is qualitatively represen- 
tative also for the case of differential rotation. (2) The 
dynamo-originated strong magnetic fields inside strange 
stars or neutron stars should be included in investigations 
of the r— mode instabilities in pulsars. Although a very 
short period P will be favourable for the occurrence of 
r— mode instability, a smaller P would also be favourable 
for the generation of a stronger magnetic field by dynamo 
action, thus exerting a stabilizing influence. As the tem- 
perature decreases, the bulk viscosity also decreases and 
the r— mode instability may occur if there are no other 
strong dissipative effects. Previously, Rezzolla et al. (2000) 
have shown that the interaction of r— mode oscillations 
with the magnetic field could be important and that the 
oscillations will be inhibited if the field is of the order 
~ 10 16 (il/fiB) G, where fl is the angular velocity of the 
star and fie is the value at which mass shedding occurs. 
Since the fields of PSSs (or PNSs) depend on the nature of 
turbulence and on the rotation frequency, further study is 
needed to see whether dynamo-generated fields can affect 
significantly the r— mode oscillation. 
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